Abstract. We have developed a new evolutionary syn- 
thesis code, which incorporates the output from chemical 
evolution models. We compare results of this new code 
with other published codes, and we apply it to the irregu- 
lar galaxy NGC 1560 using sophisticated chemical evolu- 
tion models. The code makes important contributions in 
two areas: a) the building of synthetic populations with 
time-dependent star formation rates and stellar popula- 
tions of different metallicities; b) the extension of the set 
of stellar tracks from the Geneva group by adding the 
. AGB phases for rrii/MQ > 0.8 as well as the very low 
I mass stars. Our code predicts spectra, broad band colors, 
, and Lick indices by using a spectra library, which cover 
a more complete grid of stellar parameters. The applica- 
tion of the code with the chemical models to the galaxy 
NGC 1560 constrain the star formation age for its stellar 
population around 10.0 Gy. 



(N 



> 
(N 

m 

o 

O 

o 



Key words: Galaxies: Galaxy Evolution — Stellar pop- 
ulations — individual galaxy(NGC 1560) 



A&A manuscript no. 

(will be inserted by hand later) 

Your thesaurus codes are: 

missing; you have not inserted them 



SPECTRAL: New Evolutionary Synthesis code. An 
application to the irregular galaxy NGC 1560 

Gerardo A. Vazquez^, Leticia Carigi^, & J. J. Gonzalez^ 

Institute de Astronomia, Universidad Nacional Autonoma de Mexico, A. P. 70-264, D.F. 04510, Mexico 
email: gerar, carigi, jesus@astroscu.unam.mx 

Received date; accepted date 



ASTRONOMY 

AND 
ASTROPHYSICS 



1. Introduction 

The study of the history of the evolution of galaxies in- 
cludes three important issues: spectral, dynamical, and 
chemical evolution. Since these three different parts of 
the galaxy evolution are very difficult to study in sim- 
ple models. Different authors has been focused to one 
of these parts. Although several efforts have been de- 
voted by different authors (for example Tinsley 197S and 
Vazdekis (1996)| ) to build complete spectro-chemical evo- 
lution codes. It has been difficult to incorporate the 
new releases in stellar models (structure evolution, stel- 
lar yields and stellar atmospheres) to the complexity of 
galaxy evolution models. The last reason has resulted in 
the division of the study of the spectro-chemical galaxy 
evolution in two complementary ways. While the chemi- 
cal evolution models are the independent variable in the 
complete spectro-chemical evolution, the spectral evolu- 
tion codes, which use population synthesis or evolutionary 
synthesis are left to use simple star formation histories to 
constrain the results and to be independent models. 

In this way, we find sophisticated chemical evolution 
models. But, the most part of spectral evolution models 
consider the star formation in a simple star-burst, and 
a complicated star formation history can be modeled by 
many star-bursts of different metallicities and star for- 
mation rates. Thus, instead of studying a broad range of 
galaxies using simple star-burst models, it is preferable to 
develop spectro-chemical evolution models of neighboring 
galaxies using constrains to model the systems, and after 
that use the models to galaxies at higher red-shifts. 

The new code (SPECTRAL) presented here has two 
important contributions to the evolutionary synthesis 
models: the first one is the consideration of a more com- 
plete set of evolutionary tracks from the Geneva group. 



adding tracks from Chabrier & Baraffe (1997) to the very 



low mass range; the second one is the construction of syn- 
thetic populations following the chemical evolution model 
in order to derive the spectral properties for systems with 
arbitrary enrichment history. 



In this paper we initially use simple star formation sce- 
narios to test the new evolutionary synthesis code, and 
then we proceed to more complete chemical evolution 
models to constrain the age of stellar population in the 
irregular galaxy NGC 1560. 

In Section ^ we illustrate how the tracks are assembled 
to built a complete set of tracks for either the high mass 
loss rates, or the normal mass loss rates, and we describe 
the spectra library used in our code. In Section ^, we de- 
scribe the process followed to built the evolutionary syn- 
thesis code, and the transformation to the observational 
plane of the variables. In Section ^ we test and compare 
SPECTRAL with others codes commonly used in the field 
and we present a first application by predicting spectral 
variables for the irregular galaxy NGC 1560. Finally, we 
discuss our conclusions. 



2. Stellar Tracks 

We have incorporated 
the Geneva group: 



all the tracks from 



Schaerer et al. (1993a, II) 



gchaerer et al. (1993b, IV) 



Schaller et al. (1992, I) 



Charbonnel et al. (1993, III) 
[Meynet et al. (1994, V) 



Charbonnel et al. (1996, VI), Mowlavi et al. (1998, VII) 



Send offprint requests to: Gerardo A. Vazquez 



and Charbonnel et al. (1999, VIII); together with tracks 
from Chabrier & Baraffe (1997, CB) for the very low 
mass stars, in order to assemble a complete dual grid of 
evolutionary tracks for the spectral code. One set with 
higher stellar mass loss and the other with normal mass 
loss as showed in Table ^ All these tracks are pub- 
lished in such way that each evolutionary point 
over each track has an equivalent evolutionary 
point in the rest of the tracks except for tracks from 
CB. This can be seen from the first line in the Table |^ 
because the tracks have the same number of evolutionary 
points (51). This is a very important aspect, because it is 
the basis of our evolutionary synthesis code. 

We have chosen tracks with overshooting to be assem- 
bled by using linear and bilinear logarithmic interpolation 
when it is possible and the final result is shown in Table 
||. The process to assemble the tracks is described below, 
and all parameters are logarithmic. 
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Table 1. Ingredients of the assemble of tracks. 



Tracks 


Phases 


Points 


Low 


High 


Motalli- 


Group 






Mass 


Mass 


citics 


I, II 


MS, RGB, 


51 


0.8 


120.0 


a 


III, IV 


hbS AGB^ 










V 


MS, RGB, 


51 


12.0 


120.0 


b 




HB, AGB 










VI 


MS, RGB, 


54 


0.8 


1.7 


c 




HB, AGB 










VII 


MS, RGB, 


51 


0.8 


60.0 


d 




hbS AGB^ 










VIII 


MS, RGB, 


51 


0.4 


1.0 


e 


CB 


PMS, MS, 


4 


0.075 


0.8 


/ 


a 


Z = 0.001, 0.004, 0.008, 0.02, 0.04 & M oc 


^0.5 


b Z -. 


= 0.001, 0.004, 0.008, 0.02, 0.04 & 


: M oc 5.0 






c Z = 0.001, 0.02 & 


M oc 


^0.5 






d Z 


= 0.1 & M 


ocZ°- 


5 






e Z = 0.001, 0.02 & 


M oc 


ZO.5 





/ Z = 0.0002, 0.00063, 0.001, 0.002, 0.0063, 0.02 & M = 
^Only for stars /Mg > 1.7 

Table 2. Main features from the new set of tracks. 



Tiacks 


Phases 


Points 


Low 


High 


MetalU 


Group 






Mass 


Mass 


cities 


Geneva 


PMS- 


84 


0.08 


120.0 


a 


+ Lyon 


AGB^ 










Geneva 


PMS- 


84 


0.08 


120.0 


b 


+ Lyon 


AGB^ 











a Z = 0.001, 0.004, 0.008, 0.02, 0.04, 0.1 & M oc Z° 
b Z = 0.001, 0.004, 0.008, 0.02, 0.04 h M oc 5.0 x Z' 
^PMS-AGB means PMS, MS, RGB^, HB^, AGB^ 
^Only for stars mi/Mg > 0.7 



First of all, we completed phases covered at each metal- 
licity. This means that tracks, which do not have post- 
main sequence evolution for masses in the range Wi/M© < 
1.7 {rrii is the initial mass, mass on main sequence) were 
completed to get 51 evolutionary points covering phases 
from the post-main sequence RGB, HB, and AGB. 

We have used a IDL code to visualize each interpola- 
tion in tracks, using weights for linear interpolations de- 
fined: 

^/i = 77— 7^, w/o = 1.0-u;/i, (1) 

where fx is the desired variable for which the properties 
are calculated when we interpolate in mass, metallicity or 
evolutionary phase. The sub-indices mean adjacent vari- 
ables in tracks for the interpolation. New variables {V) 
were calculated using the equation: 

Vif) = X wfo + V{f,) X wh. (2) 

Phases beyond RGB were not present for stars in the 
mass range 1.0 < rrii/MQ < 1.7, and post-main sequence 
phases were not present for stars in 0.8 < rrii/MQ < 1.0 at 
some metallicities. The first step in the new set of tracks 
was to complete tracks up to RGB for those masses in 
0.8 < nii/MQ < 1.7 using the last equations. 

The second step was to add phases as HB and AGB for 
the range in masses 0.8 < rrii/MQ < 1.7. To do the inter- 
polation and extrapolation in this phases we used tracks 
from two sets at metallicities Z — 0.001 and 0.02 men- 
tioned in the tracks group VI in Table |^. In this way, it 
was necessary to interpolate and extrapolate this phases 



in metallicity. We followed in this part the shape of tracks 
in groups I, II, HI, IV, and VII for the RGB to add the 
other two phases HB and AGB in the next form: 

We completed the evolutionary points for the RGB in 
the tracks group VI up to 30 instead of original 3 using 
the fact that luminosity grows in this phase linearly as 
explained in I. We did the same for the rest of tracks for 
masses nii/MQ > 1.7 in groups I, II, III, IV, and VII to 
get 81 evolutionary points for tracks instead of 51 in the 
first group. 

Furthermore, to interpolate and extrapolate the post- 
RGB phases for the rest groups in metallicity we use the 
next method: 

The metallicity against the shape of each stellar vari- 
able {age, current mass {rric), L/Lq, and Teff) is fitted 
for linear functions fi{Z) and f-ziZ) for both groups of 
tracks I, II, HI, IV, VII, and VI, respectively, by using the 
metallicities at which both groups coincide, in this case 
Z =0.001 and 0.02. 

Then, we defined Q{Z) = hi^) I !i{Z), as the factor 
to convert the shape of fi{Z) in to f2{Z) shape, or vice 
versa. Actually, we used the last point in RGB for the first 
tracks group as the first point in the new set of tracks and 
we interpolate in metallicity A{age), A(mc), A{L/Lq), 
A{Teff) from the second group of tracks with the last 
method. Then, the second point in the new set of tracks 
is the last point in the RGB plus the first A interpolated, 
and so on. 

The final step to assemble the tracks was to incorpo- 
rate the very low mass stars by using models from CB. 
An interpolation in metallicity has been done to put these 
tracks as the same metallicities than Geneva tracks. Then, 
we have identified over the four evolutionary points of 
tracks, ZAMS points and pre-main sequence points to ob- 
tain in the full new set of tracks 84 equivalent evolutionary 
points from pre-main sequence up to AGB. 

For comparison we show in the Figure |l| and Figure |^ 
Geneva tracks and the new set of tracks with new phases, 
and new stellar tracks for very low masses with normal 
mass loss rate, respectively. 

We can see in Figure || the pre-main sequence phases 
included for masses rrii/MQ < 0.7. In the same way, tracks 
for Z — 0.1 are very poorly interpolated because these 
tracks are originally incomplete, and have been completed 
with tracks with Z = 0.04 for masses rrii/MQ > 60. Fi- 
nally, it has been obtained two homogeneous grids to inter- 
polate the stellar variables in time, mass, and metallicity 
either, one with normal mass loss (ill cx and second 

one using tracks from the V group for masses in the range 
12.0 < uh/Mq < 120.0 in the Table with high mass loss 
(Af oc 5 X and stars in the range of rrn/MQ < 12.0 

completed in the same way as the set with normal mass 
loss showed in Figure ^. 

3. The Evolutionary Synthesis Code SPECTRAL 
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Fig. 1. Log(Te//) vs. log(L/LQ). Tracks I, II, III, IV, and 
VII from Geneva group, for masses 0.8 - 120.0 Mq. 




Fig. 2. Log(Te//) vs log(L/L0). New sets, of tracks, for 
masses 0.08 - 120.0 Mq. 



3.1. The Input and Synthetic Stellar Populations 

The evolutionary synthesis code has been developed in 
such away that it is possible to use the output of chemical 
evolution models taking from these models the foUowing 
parameters: 

-Duration of star formation, tspH- 

-Star Formation Rate, "^{t). 



-Metallicity of interstellar gas at time t, Z{t). 

The code needs the number of stars n{mi, t) formed at 
any time in the interval mi < m^/AfQ < mh defined as 
follows: 

n{m^,t)^^{m^)^{t) (3) 
where mi, m^ are the limits in mass and ^(mi) is the 
initial mass function. 

Z{t) is the metallicity of gas at which a new stellar 
population is born. Therefore, the new stars have a differ- 
ent metallicity or at least equal, predicted by the chemical 
evolution model. 

When we see the synthetic stellar population formed 
over the whole continuous star formation history from t = 
up to tsFH at Age of 10.0 Gy for instance, we need to 
calculate how many stars are still alive since they have 
formed in each step of time with the SFR (in this case 
those stars in the tip of AGB). If tspH = Age then the 
result of this procedure is that most of the stars formed in 
the first generation have died after 10.0 Gy, and just the 
very low mass stars are alive. The last generation formed 
at 10.0 Gy is complete and the age for this population 
(As) is zero. This means: 

As = Age - t (4) 
where t runs from 0.0 up to tsFH which could be the same 
as Age. 

The final result at 10.0 Gy is a mix of populations after 
the chemical evolution has regulated the SF. 

One of the problems in the evolutionary synthesis 
codes is the oscillations in the evolution of spectral pa- 
rameters, because the post main sequence phases are not 
well sampled, due to the faster evolution of stars after the 
hydrogen burning phase. As shown in the Section H, the 
problem is solved by interpolating in phases (evolutionary 
points) over the post hydrogen burning evolution, and in- 
terpolating in time and mass over main sequence phases. 

To calculate which stars are alive at Age in the post- 
main sequence region, we first calculate the masses of stars 
at each evolutionary points since the turn off point up to 
AGB using relations ta{Z,mi,j) vs (where ta is the 
age of each evolutionary point j over each track and initial 
mass mi). If those ages satisfy the condition 

ta{Z,mi,j) < As< ta{Z,m^,j + 1), (5) 
this means, 70 ages every time step, which cover exactly 
important post-main sequence evolutionary phases. The 
interpolation in phases is extended to metallicity using 
equation (1) and (2) to obtain: 
ta{Z,mi, j) — ta{ZQ,mi, j) x wZq+ 

ta{Zi,mi, j) X wZi, 
which represents the age of each 70 evolutionary points ta, 
at required metallicity for all masses in tracks m^. Masses 
at each 70 evolutionary points at required Ag are calcu- 
lated using equation (1) for the time ta(Z,mi) in the fol- 
lowing way: 

mi{Z,ta) = mi{Z,mig,ta) x wto+ 

mi{Z,mij^,ta) x wti 
where weights in time {wto y wti) use ta{Z, mi, j) and As. 



(6) 



(7) 
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Table 3. Features of SED's and broad band colors cali- 
bration used in this work. 



Para- 


Lcjcunc ct al. 


meters 


(1997, 1998) 


T.ff (K) 


(2000, 50000) K 


logs 


(-1.02, 5.5) 


[MIH] 


(-5.0, 1.0) 


AA (nm) 


(9.1, 160000) 


n(AA) 


1221 


Colors 


14 



To sample those masses in main sequence and pre-main 
sequence phases, we divide the mass range mi < nii < 
rriToff, where rriToff is the mass in the turn off point, in 
80 masses to get a 150 masses for the synthetic population 
for each time step of evolution. 

Once obtained this masses, we can calculate the stel- 
lar properties like jtLc, L/Lq, Teff interpolating in time, 
mass, and metallicity by using the equations (1) and (2) 
shown in Section 0. 




SPECTRAL01 



-3 -2 -1 

log(Age[Gyl) 



3.2. Transformation from theoretical to observational 
plane 

The stellar parameters in the synthetic population 
are transformed to the observable parameters by 
using the calibrati on of Lejeune et al. (1997) , and 



Lejeune et al. (1998) which includes low resolution spec- 



tral energy distribution's (SED's) and broad band colors 
Table ^ shows the grid of parameters covered by this li- 
brary. 

Following the same procedure as in tracks, we inter- 
polate linearly over this grid using the same equations 
described in the Section g, but now for [M/H], Teff, 
and log{g) (gravity), where [M/H] — log{Z/ZQ), with 
Zq = 0.02. Besides the spectrum, a nd broad band colors, 
we have included the Lick indices Worthey et al. (1994) 
using the routine developed by J.J. Gonzalez in 



Gonzalez (1993; 



To integrate the observable variables for all alive stars 
we have to transforms in light or flux those variables given 
in magnitudes like colors an d indices. In this way, we use 
the formalism described in Gonzalez (1993) to integrate 
the Lick indices, and 

L/,=%„xlO-«-4x(/-/o), ^ (8) 

where fi can be the V filter, and fo the bolonietric correc- 
tion, for instance, and we follow the same procedure for 
the rest of filters using the rest of colors. 

The spectral variables are weighted by the star number 
and the average values are: 

rAge jrnit) ^^^^^ ^ 



< Lf >-- 



r 

Jo 



(9) 



C'Q'>w{m',t')dm'dt' 
where the weights to obtain the average are defined 

n{m',t') X Fc For indices 
w{m',t') = ^ n{m',t') x Lboi For colors (10) 
n{m/,t') For spectra. 



Fig. 3. Predicted color index evolution of star-burst at 
Zq by three different codes; PEGASE97, GISSEL96, and 
SPECTRAL. 

and, Lf{Z,m\t') is the stellar variable to integrate, Lboi 
is the bolometric luminosity and m{t) is the mass of star, 
which is alive at time t. 

Finally, integrated variables are obtained by using : 

< (/i - fo) -2.5 X /og(^i^), (11) 

for indices and colors, and 

< magf -2.5 x log{< Lf >), (12) 
in the case of magnitude indices. 

4. Evolution Models 

4.1. Tests for the Evolutionary Synthesis Code 

We have used instantaneous star-burst models to test and 
compare our code with others commonly used. 

In the first part, we test colors and spectral indices; in 
the second part spectra for young and old populations. 

4.1.1. Colors and Spectral Indices 

For the first part we assume the next set of conditions: 

An ^(mi) oc m", with a = —2.35, low mass limit, 
mi ~ 0.1, and high mass limit, m^ = 100.0; and for SPEC- 
TRAL, we use Geneva tracks with normal mass loss rate. 

1) In the first model we compare the evolution 
of colors a,t Z = Zq with those obtained from PE- 



GASE97 (|Fioc fc Rocca-Volmerange 1997| , FR) and GIS- 
SEL96 ( [Bruzual fc Chariot 1993| , BC). The Figure | 
shows this comparison. 

We can see from the Figure |[ that there is a good 
agreement between the three models for colors. The dif- 
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Fig. 4. Colors evolution as Figure || for star-bursts at Fig. 5. Spectral indices evolution as the Figure ^ for star- 
different metallicities obtained with GISSEL96 (dashed bursts at different metallicities. GISSEL96 (dashed Unes), 



hues) and SPECTRAL (fuU lines). 



and SPECTRAL (fuU hues). 



ferences are due to the use of different spectral library 
(PEGASE97), or different stellar tracks (GISSEL96). In 
the same Figure || it is possible to see a pick in the 
(y — K) color, around log{Age[Gy]) = —1.0, which is 
produced by thermal pulses generated by stars in range 
(5.0 < TUi/MQ < 10.0), explained by FR. For a young 
population (around log{Age[Gy]) — —2.0) we can see an 
agreement between our model and that from FR in col- 
ors, and the agreement occurs with colors from BC in the 
inverse case when we compare old populations (around 
logiAge[Gy]) = 1.0). 

2) In the second model we compare the same colors 
showed in Figure || for different metallicities with those 
obtained from GISSEL96 (Figure H). In the same way we 
compare the evolution of two Lick indices iJ^ y Mgi, in 
Figure |[ 

In the plot (a) of the Figure ^ we have compared mod- 
els for the lowest metallicity reached linearly by our code 
with models at lowest metallicity in Padova tracks used 
by GISSEL96. There is a good agreement between models 
with different metallicities except for those at Z = 0.1 in 
plot (f) because the original tracks are incomplete. 

In the case of indices, our models predict a flat 
and constant evolution for young populations (around 
log{Age[Gy]) < —1.0), while indices from BC present vari- 
ations. The differences are due to they extend linearly 
to high temperatures the fitting functions (depending on 
Teff, log{M/ H), and log{g)), which predict the value of 
each index, and we use the original fitting functions with 
temperature limit around 13,000 K. There is an agree- 



ment between both models for populations in the range of 
log{Age[Gy]) > -1.0. 

4.1.2. Spectra 

We have tested spectra in a different way. We have divided 
the test for two ranges, for young stellar populations and 
for old stellar populations. Both ranges were subdivided 
in three metallicities. We use instantaneous star-bursts as 
the star formation rate. 

For young populations we assume the next set of 
conditions: 

An ^(mi) cx m", with a — —2.35, low mass limit, 
mi = 1.0, and high mass limit, nih = 100.0, and for SPEC- 
TRAL we use Geneva tracks with high mass loss rate. 

1) In the first model we compare our spectra with 
different ages in the range of 10^"^ < t/Gy < 0.9 
at Z = Zq with spe ctra produced by STARBURST99 
( iLeitherer et al. 1999) ). Figure shows quantitative differ- 
ences taking the log^Lsgg^/Ls,^) + 0.5k {k = 1, 13), where 
L599 is the luminosity from STARBURST99's spectra and 
Ls from SPECTRAL. AU models from STARBURST99 
used here are without nebular emission. 

There is a very good agreement, and it is comprehen- 
sive, because both of codes use the same spectra library 
and the same stellar tracks. There is a contribution in the 
UV part of spectra for ages at 3 and 5 My because the WR 
phases at that ages, that is the reason why the luminosity 
increase in UV part and the rate of luminosities between 
both sets of models is higher. The same features are not 
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LOG(Wavelength [Angstroms]) 



Fig. 6. Quantitative differences between our spectra and 
those from STARBURST99 at different ages, for a star- 
burst at Z = Zq. 



6 



4 



2 











1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


- 


\ \ 
\\\ 








\ 

A 


^^^^^'''"^^ """^ 200My 


-""""'^ 100My 


' 50 My 




30M/ 




20My 




\_\ 


10My 


r 








/^ii^^V-p^ 2My ^^--""^^ 
TMy~ — — — „„ . 







2.0 2.5 3.0 3.5 4,0 4.5 5.0 

LOG(Wavelength [Angstroms]) 



Fig. 8. Quantitative differences between our spectra and 
those from STARBURST99 at different ages, for a star- 
burst at Z = 0.04. 




2.5 3.0 3.5 4.0 4.5 

LOG(Wavelength [Angstroms]) 



5.0 



Fig. 7. Quantitative differences between our spectra and 
those from STARBURST99 at different ages, for a star- 
burst at Z = 0.001. 



visuahzed in our spectra because the hbrary for massive 
stars with strong stellar winds arc not considered in our 
spectra library. This increase of luminosity in UV part ap- 
pears for STARBURST99 models at high metallicity, but 
not for low metallicity and the rate of luminosities is lower 
for low Z . 

2) In the same way, the second model compares spectra 
with different ages at Z = 0.001. Figure shows quanti- 
tative differences of these models. 



3) The final comparison is done for spectra a,t Z = 
0.04. Figure ^ shows quantitative differences. 

These three sets of models show a very good agreement 
with models from Leitherer et al. (1999). The same fea- 
tures for the case at Zq remain in models (2) and (3). This 
agreement shown in three different star-burst populations 
at different metallicities is present in the rest of metallic- 



ities of Geneva tracks as showed in Vazquez (2001 



For old populations we assume the next set of con- 
ditions: 

An ^(nii) oc m", with a = —2.35, low mass limit, 
mi = 0.1, and high mass limit, rrih = 100.0; and for SPEC- 
TRAL we use Geneva tracks with normal mass loss rate. 

1) In the first model we compare our spectra with dif- 
ferent ages in the range of 1.0 < t/Gy < 20.0 a.t Z = Zq 
with spectra produced by GISSEL96 (BC). Figure ^ shows 
quantitative differences taking the log(La96k/ Lsk) + k 
{k = 1,6), where Lqqq is the luminosity from GISSEL96's 
spectra and Ls from SPECTRAL. 

It is easy to see some differences from Figure ^, which 
are due to the stellar tracks used for each code, Padova 
(GISSEL96) and Geneva (SPECTRAL), respectively. The 
considered spectral library is the same, but not the same 
version, which is a source of differences as well. The main 
differences are in the blue part of the spectra, where we 
find a contribution in UV due to old stars in phases like 
HB and planetary nebula predicted by additional stellar 
library used by GISSEL96, and there is not a counterpart 
in our models. Again, this is the reason why the rate of 
luminosities is higher for those models from GISSEL96. 

2) In the same way, the second model compares spec- 
tra with different ages at ^ = 0.004. Figure |l^ shows 
quantitative differences of these models. 
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Fig. 9. Quantitative differences between our spectra and Fig. 11. Quantitative differences between our spectra and 
those from GISSEL96 at different ages, for a star-burst at those from GISSEL96 at different ages, for a star-burst at 
Z = Zp,. Z = 0.05. 
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Fig. 10. Quantitative differences between our spectra and 
those from GISSEL96 at different ages, for a star-burst at 
Z = 0.004. 



3) The final comparison is done for spectra at Z — 
0.05. Figure [Tl| shows quantitative differences. 

The features shown in models (2) and (3) are explained 
in the same form as those with Zq . The rest of models at 
different metallicities are shown in [Vazquez (200l") . 



4-2. Models for the irregular galaxy NGC 1560 

Carigi, Colin & Peimbert (1999) discussed the chemi- 
cal evolution of three irregular galaxies: NGC 1560, II 



Zw 33 and a mean irregular galaxy. In the attempt to 
apply spectro-chemical evolution models to this three 
galaxies, we look for spectro-photometric measurements 
of these galaxies in the literature. In this work, we com- 
pute spectro-chemical evolution models only to NGC 1560 
because: i) For NGC 1560 and II Zw 33 the nonbaryonic 
dark matter inside the Holmerg Radius is known, which 
constrains the chemical evolution models, ii) For NGC 
1560, but nor for II Zw 33, neither the sample galaxies 
that form the mean irregular galaxy, there are photomet- 
ric measurements to constrain the spectral models. 

We have used close box models 

Peimbert (1999)| 



from 
dated 



Carigi, Colin 



up- 



with yields from Maeder (1992) and 



van den Hoek & Groenewegen (1997) to be consistent 



with stellar evolution models with normal mass loss rate 
used by our code. Chemical evolution models reproduce 



two constraints O/H abundance (Richer & McCall, 1995) 
and ^ = Mgas/Mtotai ( [Walter et al. 1997| [Broelis 1992] ) 



Three models with different ages and star formation rates 
to match the gas metallicity observed are not able to give 
us information about the age of star formation history. 
Colors of this galaxy may constrain the age if we use 
evolutionary synthesis. 

These three chemical models results were used by our 
code to constrain the stellar population age which domi- 
nates the light in the irregular galaxy NGC 1560. 

These models use the basic form for the IMF from 
Kroupa, Tout & Gilmore (1993), in the range of 0.01 < 
m/MQ < 120.0 as follow: 

if 0.01 < mjMQ < 0.5, 
if 0.5 < rui/MQ < 1.0, (13) 
if 1.0 < rrii/AlQ < 120.0 



$(TO.i 
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Fig. 12. Broad band colors evolution for the three models 
at 0.1 Gy (dotted line), 1.0 Gy (dashed line), and 10.0 Gy 
(continuous line). 



where a depends on chemical constraints and age of the 
model, and had three different values for every model: 

( -2.21 0.1 Gy Model 
a = < -2.24 1.0 Gy Model (14) 
[ -2.27 10.0 Gy Model 

We use $(mi), ^'(t), and Z{t) from chemical models 
in SPECTRAL to produce a set of models. The set was 
produced by using tracks with the normal mass loss rate. 
Figure |l2| shows the evolution of some broad band colors 
predicted for this galaxy. 

The rest of broad band colors and Lick index for this 
set of models are shown in 



Vazquez (2001), The mean 



spectra produced for the models are shown in Figure |13|. 

The main goal in this work is to compare our predic- 
tions for this galaxy with those obtained observationally. 
For this galaxy we have obtained {U — B)o and {B — V)o 
from the |de Vaucouleurs RC3 , which are mean values. The 
evolution and the comparison is shown in Figure |l^. 

Figure |l^ shows that the model which has a star for- 
mation age of 10.0 Gy fits in the best way the observations, 
it does not occur with younger models. 

A better fit would be obtained if there were more col- 
ors, indices or the full spectrum in the literature. 

5. Conclusions 

We have extended the tracks from the Geneva group, pro- 
ducing two sets of tracks, one with high mass loss rate 
for masses in the range of 12.0 < rrii/MQ < 120.0 and 
the other with normal mass loss rate in the range of 
0.08 < rrii/MQ < 120.0. These sets of tracks are complete 



> 
I 

oa 



10.0 Gy 




1.0 Gy 



0.1 Gy 



(U-B) 

Fig. 14. {U — B) vs. {B — V) evolution for the three mod- 
els are compared with observed values. Dotted line shows 
the 0.1 Gy model, dashed line for the 1.0 Gy model, and 
continuous line for the 10.0 Gy model. Box shows the error 
for the values from the literature. 



and consistent for metallicities in the range 0.001 < Z < 
0.04, they cover phases from main sequence up to asymp- 
totic giant branch in the range 0.8 < rrii/MQ < 120.0, 
and pre-main sequence and main sequence for masses in 
0.08 < rrii/MQ < 0.7. Together with the spectral library, 
the code SPECTRAL is able to predict the spectral prop- 
erties from stellar populations under a wide variety of con- 
ditions. 
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1.0 Gy 



A 
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Fig. 15. [U — B) vs. {B — V) from models compared with 
values from the literature. 



With this code, it is possible to build synthetic pop- 
ulations as the sophisticated chemical evolution models 
predict, as long as metallicity is in the range 0.001 < Z < 
0.04. 

Results for star-burst evolution in broad band colors, 
spectral indices and spectra from our code fit very well 
with the same results obtained with other codes. SPEC- 
TRAL could be used in young stellar populations like in 
star-burst or blue galaxies, or stellar populations inside 
cores of AGN's. In the same way, the code could be used 
for old stellar populations like those in early type galaxies 
and all galaxy types in the middle of young and old stellar 
populations. 

In applying just chemical models using observational 
constraints mentioned here for NGC 1560 it is not easy 
to distinguish what it is the age of the population. But 
using observed colors for this galaxy a spectro-chemical 
evolution model can be fitted and a population around 10 
Gy and mean metallicity of stars Z — 0.002 are predicted 
with our model. 

In the same way, from Figure the evolution of colors 
for the 10.0 Gy model is possible to see that colors match 
the box error at 5 Gy. Furthermore, we might say that 
star formation history age for the mean stellar population 
is between [5.0,10.0 Gy], though we do not have chemical 
model for 5 Gy. 

We have obtained success in reproducing the values for 
many simple star formation scenarios, and finally, we have 
matched very closely the values observed for two broad 
band colors in NGC 1560. 
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